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ABSTRACT 
The  effects  of  heat  transfer  on  the  kinetic  energy  tendencies  of  a 
perturbation  superimposed  on  a  basic  zonal  current  are  investigated.   The 
rate  at  which  kinetic  energy  is  gained  in  both  the  adiabatic  and  diabatic 
models  increases  with  time  for  medium  wave  lengths,  but  the  rate  of  in- 
crease is  less  for  the  diabatic  case.   At  the  longer  wave  lengths  the 
energy  change  for  the  diabatic  model  shows  a  decrease  with  time  becoming 
negative  toward  the  end  of  the  time  interval  investigated.   The  adiabatic 
model,  in  contrast,  has  zero  energy  tendency  at  long  wave  lengths.   The 
rate  of  kinetic  energy  change  increases  as  larger  values  of  thermal  wind 
are  selected,  and  for  the  adiabatic  model  the  wave  length  range  for  which 
energy  changes  occur  also  expands  with  increasing  thermal  wind.   The 
energy  gain  initially  is  a  maximum  when  the  thermal  wave  lags  the  pressure 
wave  by  a  phase  angle  of  90  degrees  and  the  maximum  tendency  to  lose 
energy  occurs  when  the  lag  is  270  degrees. 
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Jnt   =  the  heat  gain  or  loss 
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©C  =  the  relative  phase  angle  between  temperature  and  pressure  wave 
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I.   Introduction 

To  date  simple  baroclinic  models  have  been  unable  to  improve  signif- 
icantly on  the  results  obtained  by  one  parameter  models  in  the  prognosis 
of  mid-tropospheric  circulation.   Recent  investigations  have  shown  that 
the  incorporation  of  previously  neglected  terms,  which  were  felt  to  be 
negligible  for  short  range  forecasting,  can  significantly  alter  the 
characteristics  of  synoptic  waves  produced  by  baroclinic  forecasting 
models,  thus  bringing  these  waves  into  closer  proximity  with  snyoptic 
waves  observed  in  nature. 

•In  one  of  these  investigations,  Haltiner  f2 j  considered  the  influ- 
ence of  sensible  heat  exchange  on  the  properties  of  harmonic  perturba- 
tions superimposed  on  a  zonal  current.   The  study  revealed  that  this 
type  of  heating  reduces  the  instability  of  short  and  medium  waves,  but 
that  the  instability  of  long  waves  ia  increased. 

The  purpose  of  this  work  is  to  compare  the  kinetic  energy  changes 
averaged  over  the  interval  of  one  wave  length  in  a  baroclinic  model  with 
heat  transfer  present  to  one  when  heating  is  absent.   In  particular,  the 
manner  in  which  heating  alters  the  instability  of  long  waves  will  be 
examined.   A  simple  two-level  quasi-geostrophic  model  is  constructed   * 
and  an  equation  expressing  the  time  rate  of  change  of  kinetic  energy 
in  a  harmonic  perturbation  superimposed  on  a  zonal  current  is  derived. 
The  kinetic  energy  is  computed  numerically  for  the  range  of  wave  lengths 
corresponding  to  synoptically  observed  waves.   At  each  wave  length  com- 
putations are  made  for  varying  values  of  thermal  wind  and  phase  lag  be- 
tween the  pressure  wave  and  the  temperature  wave.   The  change  of  kinetic 
energy  with  time  is  evaluated  in  time  steps  of  three  hours  over  a  thirty 
hour  period  and  the  value  of  the  tendency  for  the  heating  and  the  non- 


heating  model  are  compared. 

In  addition,  computations  are  made  of  the  values  of  C./L,  which  is 
proportional  to  the  growth  rate  of  the  wave,  and  the  values  of  phase 
speed  to  verify  that  the  results  obtained  from  the  simple  model  used  are 
in  qualitative  agreement  with  those  obtained  by  Haltiner. 

Interest  throughout  the  study  is  confined  to  unstable  waves  only 
and  solutions  which  lead  to  neutral  waves  are  eliminated. 


II.   Mathematical  Development 

The  atmosphere    is   represented   by  a  four-layer  baroclinic  model 
using  the  750   and   250  millibar   levels   as   data   levels   and   centering   atten- 
tion  at    500  millibars.      Since  the  effects   of    latitude  variation  are  well 
known,    latitude    is   held  fixed   at   45  N.      The  vertical   velocity   is   assumed 
to   be   zero   at   both  the  top  and  the  bottom  of  the  model    atmosphere,   the 
effects  of   friction  at  the  earth's   surface   are   neglected,    and  the   static 
stability  of  the   atmosphere   is   taken  to  beQ  constant. 

The  vorticit'y  equation   is    applied   at   the    levels  corresponding  to 
250   and   750  millibars    in  the   following   form: 

cJtL 

The  continuity  equation  is  used  to  substitute  for  the  divergence  in  the 
above  equation  and  then  the  derivative  of  the  vertical  velocity  with 
respect  to  pressure  is  approximated  by  a  finite  difference  ratio.   Geo- 
strophic  approximations  are  substituted  for  the  vorticity  and  wind.   The 
sum  and  difference  of  this  equation  at  the  data  levels  are  computed  with 
the  following  results: 


A 

u 


(2) 


vl(fr^n^ti^J^^, 


At  the  middle  level  (500  mb)  the  First  Law  of  Thermodynamics  is 
applied  in  the  following  form: 

9t  r 


where 


K   m  (  >P~    cf)7    . 


The  temperature  at   this   level   is   approximated   using  the  equation  of 
state  and  the  hydrostatic  equation  to   give 

This   approximation  is   used   in  O)   and  used  to   replace  Co    in   (3). 
With  the   following   substitutions: 


I 


(5) 


our  system  of  equations  becomes 

Let  the  basic  flow  and  the  surface  be  considered  to  be  in  thermal 

equilibrium.   Similar  to  Doos  £lj  ,  the  exchange  of  sensible  heat  is 

then  assumed  to  be  proportional  to  the  difference  between  the  surface 

perturbation  air  temperature  (T  )  and  the  temperature  of  the  ground  (T  ) 

a  g 

and  a  decreasing  power  of  pressure  in  the  following  manner: 

where  G  is  a  proportionality  constant. 

Since  ridges  in  the  thermal  wave  will  correspond  to  heat  sinks  for 
the  air,  and  troughs  will  be  heat  sources,  let  the  transfer  of  heat  be 
proportional  and  opposite  in  sign  to  the  perturbation  thickness.   There- 
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Kir  -  -<vi.  5J 


fore 

°r^  dt  (8) 

where  K  is  a  proportionality  constant  with  units  sec"1. 

To  linearize  the  system  (5)  and  (6),  assume  that  the  solutions  are 
of  the  form 


ikLn-ct)  lkh-c£) 


and  that 


Dropping  terms    involving   products  of  perturbation  quantities    leads 
to   a  homogeneous   linear  system  MN=0,   where  M  is  the  matrix 


and  N   is   a  column  matrix  of     y0    and    H0  .      In  order  to   have  non-trivial 
solutions  to   this   system  of  equations  the  determinant  of  the  coefficient 
matrix  M  must  vanish.      Setting  this   determinant  equal   to   zero  yields  the 
frequency  equation 

where        __ 

3-j  at  r    ^    .    ■ 

It  can  be  seen  from  this  equation  that  the  presence  of  heating  in- 
sures that  the  phase  velocity  will  be  a  complex  number  regardless  of 
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the  wave  length.   If  the  heating  is  taken  to  be  zero,  it  can  be  shown 
that  the  phase  velocity  will  be  complex  only  for  a  certain  range  of  wave 
lengths.   The  width  of  this  range  will  depend  on  the  strength  of  the 
thermal  wind  and  the  value  of  the  static  stability. 

Using  the  primary  roots  of  the  complex  number  under  the  radical  sign, 
the  phase  velocity  has  two  solutions:  C  when  using  the  positive  sign  in 
front  of  the  radical,  and  C_  when  using  the  negative  sign.   The  complete 
solutions  for  w    and  77"  are  therefore  given  by 


h  t 


JL  4-      %     A-  (ID 


and 


where  the  primes   are  omitted  and  henceforth    lp   an^    77"* will   be   perturba- 
tion quantities  unless  otherwise  specified. 

Next  take  as   initial   conditions   at  time  t=0 

i         A     tk/  „      */&>*) 

where  A  and  B  are  real  numbers  and  °^  is  a  relative  phase  difference 
which  is  positive  when  the  thermal  field  lags  the  stream  field.   From 
(11)  and  (12)  it  follows  immediately  that  at  t=0 


tf+K**  */**;=& 


ioc  (13) 


Along  with  (13),  use  of  the  linearized  form  of  (5),  once  for  C+  and 
once  for  C_,  yields  four  equations  in  %  ,  ^  >  ^T  >  and  77"  •      Letting 
CR  equal  the  Rossby  wave  speed,  the  following  solutions  result: 
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i>f  ~    A - (14a) 

~o   -    B*   -    ^ : 


r  ft''  It"'  K    and  "i 


The  above  solutions  for  (L+ ,    lL~  •>    Jf     and  ft     are  complex  func- 
tions.  Therefore,  they  can  be  separated  into  their  real  and  imaginary 
parts 


and  similarly  for  CP  ,  7?^  and  7T  • 


As  mentioned  previously,  the  phase  velocity  can  be  a  complex  num- 
ber depending  on  whether  heating  is  present  and  the  values  of  the  thermal 
wind  and  static  stability.   The  phase  velocity  can  then  be  represented  by 


*+=  «,^<W«***<0  e  .  *,+*:*;-  («i"J) 


(15) 


where 


and 
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The  kinetic  energy  of  the  perturbation  may  be  expressed  in  terms  of 
the  square  of  the  velocity  of  the  perturbation  wind.   Since  V    is  inde- 
pendent of  y  and  considering  a  unit  mass,  the  kinetic  energy  per  unit 
mass  is 


rt-i'*-i(UT 


A  solution  for  V  is  obtained  by  substituting  into  (11)  the  solutions 
for  C+  and  C_  as  given  in  (15)  and  then  differentiating  (11)  with  respect 
to  X.   The  real  part  of  V  is  squared,  multiplied  by  one  half  (1/2)  and 
integrated  over  an  interval  of  one  wave  length  to  give  the  average 
kinetic  energy  of  the  perturbation  over  one  wave  length  as  follows: 


(17) 


where  yf-Jk.<ll't        and  77"  i-n  this  and  subsequent  equations  is  77"  =  3.14159. 

The  time  rate  of  change  of  this  kinetic  energy  is  then  obtained  by 
differentiating  (17)  with  respect  to  time  with  the  following  result  :: 


iff*  -.^fu)^::^)^^^)] 
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This  energy  change  equation  when  evaluated  at  some  time  t=t.  repre- 
sents the  tendency  of  the  perturbation  to  gain  or  lose  kinetic  energy. 
This  gain  or  loss  then  shows  whether  the  wave  is  amplifying  or  dampening. 
Since  only  unstable  waves  are  studied,  a  solution  leading  to  zero  corres- 
ponds to  a  neutral  wave. 

Haltiner  £3 J  has  shown  that  the  mean  energy  changes  over  one  wave- 
length for  adiabatic  conditions  can  be  expressed  as 


7T    ~  -/^(+.-*.)J* 


and  (19) 


U*>  „   +U(t-i)0 


at  J    l 

where  KE  and  PE  are  the  mean  kinetic  and  potential  energies.   The  quantity 
IT  "  tj)  *-s  e<lual  to  "^T  >  so  after  substituting  the  values  for  C+  and  C 
into  (12)  and  solving  for  Ol>   from  (4),  the  energy  transformation  func- 
tions on  the  right  hand  side  (19)  may  be  determined.   The  kinetic  energy 
changes  that  have  been  examined  in  this  paper  are  thus  a  consequence  of 
this  energy  transformation  from  potential  to  kinetic  energy. 


15 


III.      Results 

It    is  obvious   from  equation   (18)   that   a  determination  of  precisely 
how  heating  affects  the  stability  characteristics  of  harmonic  waves  would 
be  difficult.      Examination  of   the  frequency  equation   (10)   does   show  that 
with  heating  present   an   instability  term    (term  b)    is    always   present,    and 
also  that   heating   leads   to    an   instability  term  evolving  from  the   radical 
part   of  the  equation   (term  d) .      Term  b   always    leads   to   a  dampening  effect 
while  term  d  may  have  either  dampening  or   amplifying   characteristics. 
Therefore,   the  combined  effects  of  these  two   terms  will   determine  the 
energy  change  characteristics   of  the  heating  model.      In  the  adiabatic 
model   the  energy  change   is   dependent  on  only  term  d,   which  takes   on  a 
different  value   in  general  when  heating   is   present. 

To   arrive  at  quantitative   results   for  energy  changes    in  the  two   sit- 
uations,   representative  values   of  the  different   parameters    involved   are 
selected   and  equation   (18)    is   solved   numerically.      The    latitude    is   fixed 
at   45  degrees   north   latitude,    the  basic   flow  zonal  wind  remains   constant 
at   20  m/sec,    and  the  static   stability  also   remains  constant  with       /d-p      = 
-6K/100mb   and   so M     =   2.4  x   10"16   cm-2.      The   initial   amplitude  of   the 

stream  function  and   the  thickness   function  are  taken  to   be  equal  with  a 

10         2 
value  of   20  x   10        cm  /sec  which  corresponds   to   a   z-amplitude  of   20km. 

Since  both  the  thermal  wind   and   phase   lag   are  of   significant    importance 

in  the  energy  change,    four  values  of   thermal  wind    (5,    10,    15,    and   20 

m/sec.)    and    12   phase   lags,    30   degrees   apart,    are  used. 

It   remains  only  to   select   a  representative  value  for  the   heating 

proportionality  constant.      The  value   for  G  in    (7)    is   chosen  to  be 


£j=     IOO  ek.6,%     -J"'    •K"1  ^c"1 
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and  the  value  for  r  is  r=2.   This  is  similar  to  values  used  by  Doos  f 1  j 
and  corresponds  to  a  total  heat  transfer  of  63  cal  per  square  centimeter 
per  day  to  the  entire  column  of  air  for  every  one  degree  temperature  dif- 
ference between  the  ground  and  the  air  at  the  surface-   To  determine  a 
value  for  K,  the  initial  perturbation  thickness  is  translated  into  a 
temperature  perturbation  by  use  of  the  hypsometric  equation.   Since  the 
surface  air  and  the  ground  are  assumed  to  be  initially  in  thermal  equil- 
ibrium, this  temperature  perturbation  also  equals  the  maximum  difference 

10   2 
between  the  air  and  the  ground.   A  perturbation  thickness  of  20  x  10   cm  / 

sec  which  is  the  selected  initial  value  translates  into  a  temperature  of 

6.54K  and  corresponds  to  a  heat  transfer  of  654  ergs  g~  sec   .   A  value 

for  K  can  now  be  calculated  from  (8).   The  value  for  K  at  500  millibars 

is  then  arrived  at  by  using  r=2  and  turns  out  to  be 

K  =  2-2  <°1   X/o"^  sec'' 

Using  the  above  values,  the  first  items  determined  were  the  stabil- 
ity characteristics  of  the  waves.   As  mentioned  in  the  introduction  this 

is  the  value  of  the  complex  part  of  the  phase  velocity  multiplied  by  k, 

/  snr 

where  A  -  "JT  •   The  results  are  summarized  in  figures  1  through  4.   Figure 

1  is  for  no  heating  and  a  thermal  wind  of  lOm/sec.   Figure  2  is  again  for 
no  heating  and  now  the  thermal  wind  is  20m/sec.   Figures  3  and  4  are  the 
corresponding  situations  for  the  heating  model. 

To  gain  sufficient  resolution  between  the  values  for  each  model  at 
interested  wave  lengths  the  graphs  were  scaled  such  that  the  values  at 
short  wave  lengths  often  exceeded  the  bounds  of  the  graph.   Where  the 
bounds  have  been  exceeded  such  that  the  graph  plot  is  a  straight  line 
parallel  to  the  x-axis,  the  actual  values  at  these  wave  lengths  are 
greater  than  those  depicted  on  the  graph.   In  cases  where  the  values 
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are  significant  to  the  problem,  special  mention  will  be  made  of  the  values. 
It  should  also  be  mentioned  that  the  data  points  are  obtained  at  intervals 
of  1000  kilometers  and  between  data  points  the  computer  simply  makes 
straight  line  extrapolations. 

Examination  of  figures  1  through  4  reveals  the  following  items.   Heat- 
ing leads  to  a  smaller  amplification  factor  than  no  heating  and  a  larger 
dampening  factor.   The  values  at  5000km,  which  are  off  the  graph  for  a 
thermal  wind  of  20m/sec,  also  show  that  the  no  heating  leads  to  a  larger 
amplification  and  less  dampening.   The  effect  of  thermal  wind  is  readily 
seen  to  increase  the  instability  characteristics  even  though  the  value  of 
the  b  term  in  the  heating  model  is  independent  of  thermal  wind.   Also  it 
is  important  to  note  that  in  the  adiabatic  case  the  instability  is  con- 
fined to  a  finite  wave  length  interval  which  depends  on  the  thermal  wind 
and  static  stability.   In  contrast  the  presence  of  heating  leads  to  an 
instability  factor  at  all  wave  lengths,  however  the  amplification  ap- 
proaches zero  with  increasing  wave  length.   These  results  are  in  agree- 
ment with  similar  computations  made  by  Haltiner. 

Figure  5  depicts  the  relationship  between  the  speed  of  wave  propaga- 
tion in  the  two  models.   Ignoring  the  value  at  wave  length  zero  for  the 
heating  case,  heating  leads  to  a  decrease  of  phase  velocity  of  the  wave 
at  all  wave  lengths.   This  difference  is,  however,  quite  small  at  4-6000 
km.   This  also  is  in  general  agreement  with  results  obtained  by  Haltiner. 

Computations  of  the  time  rate  of  change  of  kinetic  energy  show  some 
interesting  features.   First  of  all,  if  the  energy  change  tendency  at 
time  t=0  is  evaluated  for  several  different  values  of  phase  lag,  the  max- 
imum positive  effect  occurs  at  a  phase  lag  of  90  degrees  and  the  maximum 
negative  change  is  at  a  phase  lag  of  270  degrees.   The  pattern  is  sym- 
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metrical  about  these  two  points  for  both  models  with  the  exception  that 
for  the  adiabatic  case  the  values  are  identically  zero  at  phase  lags  of 
0  and  180  degrees  while  in  the  case  of  heating  they  are  small  but  not 
identically  zero.   A  plot  of  the  kinetic  energy  tendency  at  time  t=0  for 
several  different  phase  angles  is  shown  in  figures  6  and  7.   As  the  time 
was  increased  from  zero  the  phase  lag  of  maximum  positive  and  maximum 
negative  energy  change  tended  to  increase,  so  that  by  12  hours  several 
situations  (wave  length  and  thermal  wind  values)  showed  a  maximum  posi- 
tive value  at  120  degrees  of  lag  and  maximum  negative  at  300  degrees. 
Also,  in  all  situations  the  magnitude  of  the  positive  rate  of  change  was 
greater  than  that  of  the  negative  change.   Again  the  graphs  indicate  that 
the  presence  of  heating  leads  to  kinetic  energy  changes  at  all  wave 
lengths  while  its  absence  restricts  the  energy  change  to  a  finite  inter- 
val.  This  phenomena  remains  true  throughout  all  situations  investigated. 

Another  point  that  should  be  brought  out  at  this  time  is  that  in 
the  derivation  of  equation  (18),  both  wave  components  of  (11)  and  (12) 
were  used.   The  energy  change  equation,  therefore,  represents  the  com- 
bination of  the  damped  and  amplifying  components.   Hence,  the  values  of 
wave  length,  thermal  wind,  and  phase  lag  which  result  in  positive  energy 
change  represent  situations  when  amplification  will  take  place  and  vice- 
versa. 

When  computations  for  the  kinetic  energy  tendency  for  the  heating 
and  non-heating  case  were  compared  at  t=0,  the  values  were  identical 
for  corresponding  situations  (with  the  exception  of  phase  lag  0  and  180 
degrees  which  was  previously  mentioned).   This  seemed  quite  surprising 
at  first  since  equation  (17)  for  the  diabatic  and  adiabatic  cases  is 
quite  different.   Upon  further  reflection  it  is  apparent  that  since 
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heating  has  had  no  time  to  exert  an  influence,  the  values  would  be  the 
same.   Figures  8  and  9  show  the  tendency  at  t=0  for  several  thermal  winds. 
On  the  curves  for  a  thermal  wind  of  20m/sec,  the  values  at  wave  lengths 
5000,  6000  and  7000  kilometers  can  be  seen  to  be  equal  in  both  cases. 

Next  the  tendency  is  evaluated  at  t=6  hours.   In  both  cases  the 
values  increased  in  magnitude  but  now  an  examination  of  figures  10  and 
11  show  that  the  values  for  the  heating  case  are  smaller  than  those  for 
no  heating.   When  time  is  allowed  to  increase  further,  the  difference  in 
tendency  values  for  the  two  cases  becomes  larger.   Figures  12  through  15 
are  for  t=12  and  t=24  hours  and  show  how  the  difference  increases.   By 
24  hours  the  values  at  5000  kilometers  are  off  the  graphs  but  the  raw 
data  still  gives  a  higher  value  for  the  adiabatic  situation. 

Values  of  energy  change  were  actually  tabulated  out  to  33  hours  at 
a  phase  lag  of  90  degrees.   In  the  adiabatic  model  the  tendency  value 
increased  out  to  33  hours  in  all  cases.  When  heating  was  present  the 
value  also  increased  with  time  for  medium  wave  lengths  with  the  exception 
that  after  30  hours  for  a  wave  length  of  5000  kilometers  and  a  thermal 
wind  of  5  and  lOm/sec,  the  value  started  to  decrease.  This  suggests  that 
not  only  does  the  tendency  increase  at  a  slower  rate  when  heating  is  in- 
troduced, but  also  that  it  reaches  its  maximum  value  sooner.   At  the 
longer  wave  lengths  the  value  of  the  tendency  decreases  with  time  such 
that  by  24  hours  kinetic  energy  is  actually  being  lost  by  the  waves  with 
wave  lengths  larger  than  10,000  to  11,000  kilometers.   It  was  mentioned 
earlier  that  the  phase  lag  of  maximum  energy  change  increases  with  time 
also.   This  change  of  phase  lag  does  not,  however,  alter  the  results  that 
heating  leads  to  smaller  values  of  tendency  as  the  time  increases. 

So  far  we  have  examined  primarily  the  spectrum  of  medium  length  waves 
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to  determine  the  effects  of  heat  transfer  on  the  energy  changes.   At  the 
long  wave  lengths  the  addition  of  heat  transfer  to  the  model  leads  to 
only  a  small  rate  of  energy  change.   This,  however,  is  significant  since 
in  the  adiabatic  case  the  energy  change  is  identically  zero  at  long  wave 
lengths.   That  heating  insured  the  presence  of  an  instability  term  at 
all  wave  lengths  was  observed  in  equation  (10)  and  again  in  the  results 
of  figures  3  and  4.   These  did  not,  however,  indicate  whether  the  effect 
was  to  amplify  or  dampen  the  wave.   Examination  of  the  kinetic  energy 
tendency  due  to  heating  clearly  shows  that  the  effect  depends  on  the  time 
and  the  phase  lag.   For  example,  the  series  of  graphs  on  energy  change 
show  that  out  to  at  least  24  hours,  if  the  phase  lag  is  90  degrees,  the 
heating  leads  to  an  amplifying  wave  at  long  wave  lengths.   Beyond  30 
hours  the  heating  leads  to  a  dampening  effect. 

The  energy  tendency  at  the  very  short  wave  lengths  was  not  examined 
comprehensively  and  no  explanation  is  offered  for  the  very  large  values 
of  energy  change  observed.   It  can  be  seen  from  the  graphs,  however,  that 
besides  displaying  very  large  values,  the  kinetic  energy  tendency  may 
change  abruptly  from  a  negative  to  a  positive  value  or  vice  versa  as  the 
time  or  wave  length  varies. 
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